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Abstract. We study the impact of transverse-momentum dependent parton distributions on 
detailed features of multi-jet final states, focusing on angular jet correlations in DIS data. 



Monte-Carlo simulations of events containing multiple jets are central to a number of processes 
of interest at LHC energies, see e.g. [I], and further references in [2]. These complex events, 
characterized by multiple hard scales, are potentially sensitive to QCD initial-state radiation 
processes that depend on the large- k± tail of partonic distributions and matrix elements [21[3j|3j. 
Such effects are not included in standard shower Monte-Carlo implementing collinear jet 
evolution. In this article we report on the study [5] of multi-jet DIS final states at HERA in the 
framework of unintegrated parton distributions. Despite the much lower energy at HERA, we 
stress that, owing to the large phase space available for jet production, the HERA data may be 
just as relevant as the Tevatron for understanding the extrapolation to the LHC of initial-state 
radiation effects [2]. We refer the reader to [2] for a recent discussion of /cj_-dependent parton 
distributions. 

The ZEUS collaboration [6] has recently presented measurements of multi-jet distributions 
associated with Q 2 > 10 GeV 2 and Bjorken-x as low as x ~ 10~ 4 , shown in Fig. Q] [6] with a 




Figure 1. (left) Bjorken-x dependence and (right) azimuth dependence of di-jet distributions 
at HERA as measured by ZEUS[6j. 
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comparison to next-to-leading-order calculations [7] . The plot on the left in Fig. [T] shows the 
x-dependence of the di-jet distribution integrated over Acp < 2n/3, where Acp is the azimuthal 
separation between the two high-E^ jets. The plot on the right shows the di-jet distribution in 
Acp for different bins of x. The data are in agreement with NLO results. However, the variation 
of the predictions from order-a^ to order-a^ is very large as x and Acp decrease. The theoretical 
uncertainty at NLO thus appears to be sizeable, and likely underestimated by the dashed error 
band obtained from the conventional method of varying the renormalization/factorization scale. 
The analysis [6] also indicates that for more inclusive jet cross sections perturbative predictions 
are much more stable. 

The jet correlations in Fig. [T] of [6j depend on a number of physical processes that may affect 
the stability of predictions. Part of these concern the jet reconstruction and hadronization. In 
this regard note that the jet algorithm used in [B] is expected to have moderate hadronization 
corrections [8] and to be free of nonglobal single-logarithmic components [9J. The kinematic 
cuts on the jet transverse momenta are set to be asymmetric, so as to avoid double-logarithmic 
contributions in the minimum px [H] • Nonperturbative corrections affecting the jet distributions 
at the level of inverse powers of Q are expected to be moderate for Q 2 > 10 GeV 2 . 

Further, potentially important processes are associated with higher-order radiation. Higher- 
order contributions enhanced by soft and collinear gluon emission are produced from 
configurations where the jets are nearly back-to-back [11]. To the leading-log level these 
contributions are also taken into account in standard shower Monte-Carlos, e.g. Herwig |12j . 
However, also sizeable corrections in Fig. [T] are seen in the range of decreasing Acp, where the 
two jets are not close to back-to-back, and one has effectively three well-separated hard jets. 
The corrections increase as x decreases. 

Motivated by this, Ref. [5] investigates the possibility that large higher-order terms come not 
only from the infrared emission but also from soft-gluon exchange producing relatively energetic 
multiple jets into the final state. In this kinematics, a sizeable k±is built in the initial state, 
whose effects are not fully accounted for by the collinear approximation implemented both by 
Herwig and by fixed-order calculations. Ref. [5] takes the fcj_-dependence into account both 
in the initial-state evolution and in the hard-scattering matrix elements, based on the Monte- 
Carlo implementation [13] Cascade of high-energy factorization [T3] and the fits [15] to the 
unintegrated gluon density. 

Fig. [2] shows results for the azimuthal distribution of di-jet and three-jet cross sections [5], 
compared with the measurement [6]. The description of the measurement by Cascade is good, 
whereas one can see that the collinear based parton shower with leading order matrix elements, 
as implemented by Herwig, is not sufficient to describe the measurement in the small Acp region. 
To gain insight into the physical picture of the production process it is also worth looking at the 
angular distribution of the third jet. Fig. [3] [5J shows the distribution in the azimuthal separation 
between jet 1 and jet 3, and compares it with the result from Herwig. The k± Monte-Carlo 
Cascade gives large differences to Herwig in the region where the azimuthal separations Acp 
between the leading jets are small (left hand side of Fig. [3]). It becomes closer to Herwig as Acp 
increases (right hand side of Fig. [3]), consistently with the expectation that both Monte-Carlos 
give reasonable approximations in the back-to-back region. The behavior in Figs. [2] and [3] is 
confirmed by the study of jet multiplicities [5]. The contribution of high multiplicities is found 
to be significantly increased in the /cj_Monte-Carlo Cascade compared to Herwig. A behavior 
qualitatively similar to that of the angular correlations is obtained in [5] for the distribution in 
the transverse-momentum imbalance between the leading jets. 

In conclusion, the analysis [6] of multi-jet production suggests that, while inclusive jet cross 
sections are reliably predicted by NLO perturbation theory, jet correlations are affected by 
sizeable higher-order corrections. The corresponding theoretical uncertainties are large at NLO. 
The results of [5j support a physical picture in which the large corrections are due not only to 




Figure 2. Angular jet correlations [5] obtained by Cascade and Herwig, compared with 
ZEUS data [BJ: (left) di-jet cross section; (right) three-jet cross section. 
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Figure 3. Cross section in the azimuthal angle A(^x3 between the hardest and the 3rd jet for 
small (A(p < 2, left) and large (A<p > 2, right) azimuthal separations between the leading jets [5]. 
The k±_ Monte-Carlo results Cascade are compared with Herwig. 



contributions of infrared gluon radiation but also to the growth of the fej_in the initial state, 
enhancing effects of the x <C 1 parton cascade. While the former are summed to leading-log 
level by the Herwig shower (and by the calculation [llj to the next-to- leading- log level), the 
latter require the use of shower Monte-Carlo implementing A^-dependent parton distributions 
and matrix elements like Cascade. The main signatures of the difference between the two 
showerings are found [5] in the small-A(/) behavior of jet correlations and in the jet multiplicities. 
Although this analysis is for DIS, we note that a similar situation may be expected for multi-jet 
final states at the LHC (unlike the Tevatron) , owing to the large phase space becoming available 
at LHC center-of-mass energy. 

The present k± Monte-Carlo Cascade gives already an improved description of the 
measurements compared to standard parton showers approaches. See e.g. [I] for a discussion of 
current limitations, and directions of progress, including aspects of the theory of unintegrated 
parton distributions \16\ I17j. Note that the multi-jet observables discussed in this article give 
an especially interesting probe of the approach, as they are much more sensitive than inclusive 
cross sections, and much less model-dependent than forward-region cross sections. 
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